Effect of depreciation of the public goods in spatial public goods games 
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In this work, depreciated effect of the pubhc goods is considered in the public goods games, 
which is reahzed by rescaling the multiphcation factor r of each group as r' = r{^)^ (/? > 0). It 
is assumed that each individual enjoys the full profit of the public goods if all the players of this 
group are cooperators, otherwise, the value of the public goods is reduced to r' . It is found that 
compared with the original version (/? = 0), emergence of cooperation is remarkably promoted for 
/3 > 0, and there exit optimal values of P inducing the best cooperation. Moreover, the optimal 
plat of /3 broadens as r increases. Furthermore, effect of noise on the evolution of cooperation is 
studied, it is presented that variation of cooperator density with the noise is dependent of the value 
of P and r, and cooperation dominates over most of the range of noise at an intermediate value of 
/3 = 1.0. We study the initial distribution of the multiplication factor at /3 = 1.0 , and find that all 
the distributions can be described as Gauss distribution. 

PACS numbers: 89.65.-s, 87.23.Kg, 87.23.Ge 



I. INTRODUCTION 



Evolutionary game theory [l|-|7|is widely applied to 
study the maintenance of cooperation among the self- 
ish individuals. The prisoners' dilemma game (PDG) by 
pairwise interaction ISl-flOl and public goods game (PGG) 
by group interaction have been extensively used to 
investigate the altruistic behavior. In a classical pub- 
lic goods game, N players are selected randomly from a 
large population, and each player has two choices: co- 
operation and defection. A cooperator will contribute 
a mount cost to the public pool, while a defector con- 
tributes nothing. The accumulative contribution is mul- 
tiplied by a factor r, then redistributed equally among 
all the players. In a well-mixed population, it leads to a 
rock-scissors-paper dynamics when loner strategy is in- 
troduced [13 • Considering the limitation of the classical 
game theory, some mechanisms and theoretical supple- 
ment are proposed, such as reputation and punishment 
p^l - fisj , network reciprocity [H, [T9l-l22j , voluntary par- 
ticipation [H-Ili]. 

Diversity or the inhomogeneity has been studied in 
many works, which mainly focuses on the network topol- 
ogy diversity or individual inhomogeneity [26l - l30| . San- 
tos et al. have introduced social diversity by means of 
heterogeneous graphs, and concluded that cooperation 
is promoted by the diversity associated with the num- 
ber and the size of the public goods game, as well as 
the individual contribution to each group [22]. Consid- 
ering the profit diversity of the public goods in reality, 
group diversity is introduced in the public goods game 



in which r follows a given distribution among the popu- 
lation However, in real situations, the value of the 
public goods is not invariable, but evolves because of the 
external conditions or the reasons of themselves. In this 
work, we study the evolution of cooperation in the spa- 
tial public goods games by considering the depreciated 
effect of public goods. It is assumed that each individual 
enjoys the full advantage of the public goods if all the 
people are cooperators in a single group, otherwise, the 
value of the multiplication factor of this group is reduced 
as a function of r{^Y . No qualitative change happens 
when altering the initial distribution of the cooperators. 

The paper is organized as follows. In section II, de- 
scription of the model is proposed in detail. Numerical 
simulations and the correspondent analysis are presented 
in section III. Conclusions are drawn in the Section IV. 



II. MODELS 

Public goods game is studied on a square lattice with 
periotic boundary conditions. Each player can either 
cooperate or defect, and interacts with its nearest four 
neighbors. Here each individual only participates in a 
SINGLE group. Considering the depreciated effect of 
the public goods, the multiplication factor of the group 
centered on individual x is rescaled as 
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where r is the multiplication factor indicating the full 
profit of the public goods of each group, ricx denotes 
the number of cooperators in this group, and G (=5) 
is the group size. /3 > is a tunable parameter, and 
when /3 = 0, the system returns to the original version in 
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FIG. 1: (Color online) Cooperator density pc as a function of 
multiplication factor r for different values of (3. 



FIG. 2; (Color online) Variation of pc with /3 for different 
values of r. 



which r is invariable and same in each group. The payoff 
of player x is given by 



p / "ex 



(2) 



where Sx indicates the strategy of x, Sx=l for a coop- 
erator, and for a defector. 

After each time step, a player x will update its strategy 
by choosing a neighbor y randomly from the neighbor- 
hood. We adopt the updating rule depending on their 
total payoff difference. 



W{Sx -> Sy) 



(3) 



where k denotes the amplitude of noise level. Here we 
set K = 0.1. 



III. SIMULATION AND ANALYSIS 

Simulations are carried out for a population of iV = 
200 X 200 individuals with the synchronous updating rule. 
Initially, the two strategies of cooperation (C) and defec- 
tion (D) are randomly distributed with the equal proba- 
bility of 1/2. The key quantity for characterizing the co- 
operative behavior is the density of cooperators pc, which 
is defined as the fraction of cooperators in the whole pop- 
ulation. In the all simulations, is obtained by averag- 
ing over the last 5000 Monte Carlo (MC) time steps of 
the total 45000. Each data point results from an average 
of 50 realizations. 

Figure 1 shows the variation of pc with r for different 
values of /3. One can see that pc increases with r for 



each /3, and emergence of cooperation is remarkably pro- 
moted compared with the original version (/3=0). More- 
over, /3=1.0 induces the best promotion, which suggests 
that there exit optimal values of (3 for the evolution of 
cooperation. 

Then we study the effect of /3 on evolution of coop- 
eration in Fig. 2 for different values of r. It is shown 
that for each value of r, there exit intermediate values 
of /3 leading to the highest cooperative level, and mean- 
while, optimal plat of /? is formed which broadens as r 
increases. Moreover, we can also conclude that the initial 
distribution of r has a critical impact on the evolution of 
cooperation from Fig. 2. It is shown that 5 (/3 = 0) 
or polarized {j3 — > oo) distribution has an adverse effect 
on evolution of cooperation, which is consistent with the 
results in the previous work. 

Before we further study, it is necessary to explain the 
results gained above. It is known that group diversity 
can promote cooperation in the public goods games in 
which the multiplication factor r follows a given distri- 
31| . In our model, diversity of group forms when 
according to equation (1), and meanwhile, the 



bution 

/3 > 

diversity is not fixed, but evolves with the time step as 
cooperation evolves. The spatio-temporal diversity con- 
tributes to the spread of cooperation, and thereby better 
promotes the cooperative level than the fixed spatial di- 
versity. Because in the fixed distribution of diversity, the 
players in disadvantageous group with very small values 
of r will always give up cooperating (see reference jsij). 

Figure 3 shows the time evolution of cooperation for 
different values of /3 at r=4.5. It is found that for small /3 
(/3=0.2), cooperative clusters are formed, and sustained 
as small sizes; when /3 is moderate, large clusters are 
formed, and are spreading over all the population as t 
evolves. While for the large values of /3, cooperative clus- 
ters can't be sustained, and gradually exploited by the 
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FIG. 3: Time evolution of spatial distributions of coopera- 
tors (white) and defectors (black) at r=4.5 for (a) /3=0.2, (b) 
13—1.0, (c) /3=1.8. The initial distribution of cooperators is 
uniform distribution for each condition. 
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FIG. 4: (Color online) Cooperator density pc as a function 
of noise amplitude k for different values of /? at r=4.0, 4.5, 
5.0, 5.5 respectively. 

defectors with the time step. 

We further investigate the impact of noise on the evo- 
kition of cooperation in Fig. 4 for different values of r. 
Apparently, the variation tendency of cooperator density 
Pc with K not only depends on the value of r, but also 
has the relation with the value of /3: when r is too small 
or too large, Pc decreases with k for each /3, while for 
a moderate value of r, Pc plays nonmonotonic behavior 
at some value of (5 (e.g., see r=4.5 at (3—0.2). It should 




FIG. 5: (Color online) Initial distribution of r' at /? = 1.0 for 
different values of r. Here k=0.1. 



be also noted that cooperation almost dominates at an 
intermediate value of /3—1.0 for each r. 

Given that the initial distribution of the group diver- 
sity plays an important part on the evolution of cooper- 
ation, we study the initial distribution of multiplication 
factor r' for different values of r at /? = 1.0 (see Fig. 
5) in order to figure out why /3=1.0 has such an advan- 
tageous status in Fig. 4. It is surprisingly found that 
all the initial distribution can be described as Gauss dis- 
tribution, which implies that Gauss distribution has an 
positive effect on the evolution of cooperation. 



IV. CONCLUSIONS 

In this work, evolution of cooperation was investigated 
in the spatial public goods games by considering the ef- 
fect of depreciation of the public goods. We assumed that 
each individual enjoys the full profit of the public goods 
if all the players of this group are cooperators, otherwise, 
the value of public goods is depreciated, which is realized 
by rescaling the multiplication factor r of this group to 
r{^^^)^ . So that spatial diversity of the multiplication 
factor is formed among the population, and meanwhile, 
evolves with the time step as the cooperation evolves. 
We found that compared with the original version for 
/3 = 0, cooperative level is largely improved when /3 > 0, 
and there exist moderate values of (3 leading to a best 
cooperation. Moreover, optimal flat of /3 are formed, and 
broadens as r increases. The facilitation of cooperation 
should be attributed to the spatio-temporal diversity of 
the multiplication factor caused by the depreciated ef- 
fect of the public goods: spatial diversity stimulates the 
emergence of cooperation, whereas temporal diversity of 
r contributes to the spread of cooperation. Furthermore, 
noise impact on the evolution of cooperation was studied, 
it was proved that the variation of cooperation density 
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not only depends on the values of r, but also has a close 
relation with the values of We also noted that there 
exits an intermediate value of ^ = 1.0 inducing the best 
cooperation over most of the range of noise. In order to 
explain this phenomenon, initial distribution of multipli- 
cation factor at 13 = 1.0 is studied. It was shown that 
all the distribution can be depicted by the Gauss distri- 
bution, which suggested that Gauss distribution has a 
positive effect on the evolution of cooperation. 

Our work reiterates the importance of the diversity on 
the evolution of cooperation, and meanwhile, that Gauss 
distribution has such an advantageous status may give us 



some inspiration to study the cooperation in a different 
way. 
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